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Abstract. We compute the pole mass of the gluon in QCD from the local composite operator 
formalism at two loops in the MS renormalization scheme. For Yang-Mills theory an estimate 
of the mass at two loops is 2.13Aj^g. 
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1 Introduction. 



There has been recent interest in understanding the role the dimension two gauge invariant gluon 
mass operator plays in the vacuum structure of Yang-Mills theory and QCD. For instance, see 
[UElElllElEllZlliinilinilllJIiaiiniand references therein. As was pointed out in [HI UHl ESI CH 
the perturbative vacuum one ordinarily uses for high energy computations is not stable and it 
has been suggested that in the true vacuum, various operators condense developing non-zero 
vacuum expectation values. Whilst the operator product expansion and QCD sum rules are 

usually centred on the gauge invariant operators (G^J] an d 'tpifi, where G® u is the gluon field 
strength and ip is the quark field, recent work considered the lower dimension operator \A a2 in 
the Landau gauge, El EH EEl EEU EE21 EE] , and various generalizations of it. For instance, there 
one can construct a non-local gauge invariant dimension two operator which truncates to \ A't 2 
in the Landau gauge. An example of the role such an operator can play in the infrared structure 
has been examined in ^H] where a massive gauge invariant QCD Lagrangian was studied in 
connection with vortex solutions. As \A a2 has the dimensions of a mass operator, it has been 
the subject of investigating the issue of whether the gluon can develop a mass dynamically. 
Indeed in the early work of Curci and Ferrari, |19j . an extension of this operator, which was 
on-shell BRST invariant, was included in the usual QCD Lagrangian. More recently another 
approach, known as the local composite operator (LCO) method, has been developed which 
avoids the ad hoc inclusion of a gluon mass term [HI El EH- Instead the QCD Lagrangian is 
modified to introduce an extra scalar field, a, coupled to \A a2 and with this one can compute 
the effective potential of the scalar field. It transpires that due to the development of a non-zero 
vacuum expectation value for a the gluon gains a non-zero mass in a vacuum which does not 
correspond to the (unstable) perturbative one. By contrast, applying the same formalism to 
QED, jlUj . the perturbative vaccum is stable and whilst there is another extremum where a has 
a non-zero vacuum expectation value, it corresponds to an unstable point. 

Within this formalism one can estimate the size of an effective gluon mass at one and 
two loops. In Yang-Mills theory it is of the order of 2Aj^g and is stable to the higher order 
corrections. Whilst this is roughly consistent with other estimates of a gluon mass from a wide 
range of methods (which are succinctly summarized in Table 15 of Field's article, f2Uf) the LCO 
estimates suffer several shortcomings. One of these is that the effective gluon mass used in (Hj was 
the tree object and whilst the effective potential does have the quantum corrections no account 
of the dressings of the tree quantity were included. Further, whilst all the gluon estimates are 
of a similar range, it is not clear to what extent the same mass quantity is being measured. For 
instance, in the quark sector of QCD the quark masses are all measured and compared to the 
same benchmark, which is the running mass at the scale 2GeV. This is irrespective of whether 
the pole mass of the quark was determined or, say, the running mass at another scale prior to 
using the (four loop) quark mass anomalous dimension to run the mass to the standard reference 
scale. For the same problem for a gluon mass, the anomalous dimension of the \A°^ gluon mass 
operator in the Landau gauge is now available at four loops, [2^, extending the two, |22j, and 
three loop, j2H], results which are all in the MS scheme. Remarkably the operator anomalous 
dimension is the sum of the gluon and ghost anomalous dimensions in the Landau gauge, |23U24| . 
To complete the analysis for any future gluon mass computations one requires, for instance, the 
relation between the running and pole mass of the gluon. This was initially addressed for the 
LCO formalism in j2S] where the one loop relation between these quantities was given in the MS 
scheme where the computation extended the one loop calculation of |26| for the Curci-Ferrari 
Lagrangian (with Nf = 0) itself rather than the LCO one. Moreover, an estimate of the pole 
mass was provided, |25| . by converting the effective potential of the classical gluon mass of the 
LCO Lagrangian into a potential for the gluon pole mass. Remarkably upon extremization and 
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solution of the resultant equation, a mass estimate emerged in QCD which was independent of 
the renormalization scale. For Yang- Mills theory this corresponded to 2.10Aj^g, |25| . though 
the values for Nt = 2 and 3 were significantly lower. Given the interest in gluon masses and 
the absence of relations between the various mass quantities, it is the purpose of this article to 
extend [22] to two loops by computing the gluon pole mass in the LCO formalism in the MS 
scheme in the Landau gauge. As a by-product we will be to deduce the same quantity in the 
Curci- Ferrari Lagrangian to extend the Yang-Mills result of [2^] . Another motivation for such an 
analysis, aside from determining whether the two loop corrections significantly alter the one loop 
estimates, is to ascertain whether the two loop result using the minimization criterion of j2^j 
retains the one loop renormalization scale independence. Given the fact that by analogy, except 
for the /3-function, the MS anomalous dimensions are renormalization scheme independent at 
only one loop, it would be surprising if a two loop pole mass estimate remained renormalization 
scale independent. 

The article is organised as follows. In section 2, we review the background to the problem 
including the necessary points of the LCO formalism before discussing the construction of the 
two loop gluon pole mass in section 3. Equipped with this we perform the analysis to produce 
a mass estimate in section 4 before concluding with various estimates in section 5. 



2 LCO formalism. 



We begin by reviewing the key ingredients of the LCO formalism, we will require. First, we 
define the QCD Lagrangian in an arbitrary covariant gauge as 



L QCD 



il 



(2.1) 



where a is the gauge fixing parameter, A ^ is the gluon field, c a and c a are the ghost and anti- 
ghost fields, is the quark field and the various indices range over 1 < a < Na, 1 < I < Np and 
1 < % < Nf where Nf and Na are the dimensions of the fundamental and adjoint representations 
respectively and Nf is the number of quarks. To construct the LCO Lagrangian from LQ CD we 
introduce the path integral W[J] defined by, [Hj, 



-W[J] 



VA^Dip Vtp 'Dc 'Dc exp 



d d x (l 



1 t Aa2 



U T 2 



(2.2) 



where J is the source coupled to the local composite operator \A a ^ in the Landau gauge and 
the subscript denotes bare quantities. To retain renormalizability of the action including the 
source as well as a homogeneous renormalization group equation an additional term has been 
introduced. For instance, the term quadratic in J is necessary since the vacuum energy is 
divergent as can easily be seen by power counting. This term is coupled in via a parameter 
£ and its associated counterterm <5£ is included when the action is converted to renormalized 
parameters giving 



-W[J] 



VA^VipVipVcDc exp 



d x [L 



(2.3) 



where Z m is the gluon mass renormalization constant. It transpires, [HJ, that in the LCO 
formalism one can compute the explicit form of £ by ensuring that W[J] does indeed satisfy 
a homogeneous renormalization group equation. Consequently since the coupling constant, g, 
runs in such an equation £ is constrained to satisfy a differential equation dependent on the 
/5-function and anomalous dimension of the gluon mass operator, \A a2 This equation can be 
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solved in a coupling constant expansion. As we will require the explicit form here we note that 
in QCD we have, jaUUJ, 
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(2.4) 



where TV [r a T b ^j = T F 5 ab , T a is the group generator, C A and C F are the usual colour group 
Casimirs and C,{z) is the Riemann zeta function. Consequently one uses a Hubbard-Stratanovich 
transformation to rewrite the exponential in the path integral of W[J]. This introduces the 
additional scalar field a and gives a generating functional where the source J now couples 
linearly to a field as opposed to a composite operator. Thus, |Hj, 



-W[J] 



VAJJipVipVcDcDa exp 
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(2.5) 



where L a is the LCO Lagrangian and is given by, |Hj, 



a 2 



+ 



2a 

5" aA a u A a " 



7 2 



2g 2 £(g)Zt *gc,yg)^z 

where the first few terms of £(<?) are given by (|2.4j) . 



8H{g)Z L 



(A«A a > 



(2.6) 
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Using (j2.fi |) the effective potential V(a) was constructed by standard methods at two loops in 
the MS scheme, |8"ll9*lll0|. The divergences in the vacuum energy are removed by straightforward 
renormalization since the Lagrangian, L CT , retains the renormalizability property. As the two 
loop effective potential will be required later we recall its explicit form is 
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where 
Ai = 



[13C A - 8T F N f 



-1 
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and n is the renormalization scale which incorporates the usual factor of 47re~ 7 into the MS 
renormalization scale where 7 is the Euler-Mascheroni constant. Minimizing V(a) with respect 
to the quantity a one discovers that the classical perturbative vacuum at (a) = is unstable 
and that there is a stable vacuum for a value of (a) 7^ 0. Estimates for the value of {a) were 
given in [H] by assuming that 

dV(a) 



da 







(2.9) 



5 



and then choosing the renormalization scale [i to be such that there were no logarithms in the 
equation relating the coupling constant to the value of (a). Using the renormalization group 
properties relating the coupling constant to the scale \x and hence the fundamental scale Aqqq 
in the MS scheme, the estimate for (a) was obtained which was relatively stable to two loop 
corrections. In [SJ a subsequent estimate was deduced for an effective gluon mass. This will be 
illustrated in more detail in a later section. 

However, as indicated earlier this was essentially the classical or bare gluon mass in the 
region of the stable minimum of the effective potential for a. A more appropriate quantity to 
examine would be a gluon mass where some account of the quantum corrections were included. 
In this article we will use the pole mass of the gluon as constructed from the LCO Lagrangian 
at two loops building on the previous one loop analysis of [25] . 



3 Two loop gluon pole mass. 

In this section we construct the relation between the running gluon mass of the LCO Lagrangian 
and the pole mass which is defined to be the pole of the one particle irreducible gluon polarization 
tensor. In [27] the simpler two dimensional Gross-Neveu model was studied and the relation 
between the analogous quantities was determined. The final mass estimates compared favourably 
with the known exact mass gap. Whilst we will use [2Zj as a basis for the QCD computation 
there are significant differences aside from the space-time dimensionality. The first is that L a 
has more interactions and basic fields as well as the gauge property. Second, and partly as a 
consequence of the previous point, it is not possible to fully construct the gluon 2-point function 
for all momenta and then deduce the pole of the propagator. This is also due to the fact that 
not all relevant basic 2-point two loop Feynman diagrams can be written in terms of closed 
known analytic functions for all values of the momenta. To circumvent these difficulties we 
have followed the strategy and algorithm of a similar model in the context of the weak sector of 
the full standard model. Our approach is based on the series of articles [2H1 EH1 HOI EH] which 
applies the On-Shell algorithm to the relation of the vector gauge boson poles masses in MS 
to their bare values. This package, [2E1 123!) ^ s designed to determine the value of two loop 
Feynman diagrams with massless propagators in addition to a propagator with a mass which 
is the on-shell value whose pole mass one is interested in. It uses dimensional regularization in 
d = 4 — 2e dimensions. One can extend the approach of |28[ I29| to integrals with more than 
one scale by expanding in an appropriate ratio of masses which is assumed to be small. In our 
case this complication does not occur. 

The On-Shell package, |28[ I29|. is written in the symbolic manipulation language Form, 
|32| , and for L a we have generated the relevant one and two loop one particle irreducible Feynman 
diagrams using the Qgraf package, [231 • This is converted into a Form readable format before 
applying the On- Shell procedure to determine the value of each individual diagram when the 
external momentum is set to its on-shell value. For the LCO Lagrangian we are interested in 
there are 5 one loop diagrams and 39 two loop diagrams to evaluate. 

The remaining issue is to construct the pole mass itself from the integral contributing to the 
gluon 2-point polarization. If we define the transverse part of the correction to the polarization 
tensor by 



LV(p) = U(p 2 ,m 2 ) 



(3-D 

where p is the external momentum then the pole mass is defined to be that value of p 2 which is 
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the solution to, |3*Uj . 

p 2 - m 2 - U(p 2 ,m 2 ) = . (3.2) 
If we write the perturbative expansion of the transverse part of the polarization tensor as 



(3.3) 



n=l 



then to two loops one can solve the pole mass condition iteratively to obtain the pole mass, s p , 
as, (30]) 



where 
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(3.5) 



and here m 2 = m 2 (/i) is the running mass. The actual values of Ili(m 2 ,m 2 ) are obtained 
from the On-Shell package, |281 129j . In determining the two loop part of Q3.4JI from the one 
loop diagrams, we have expanded the bare coupling constant and bare mass in terms of the 
renormalized variables before applying the one and two loop On-Shell routines. As a check 
that the final expression we obtain for the pole mass is correct we note that first the full 2- 
point function LI(m 2 , m 2 ) itself has to be finite at two loops after renormalization with the usual 
Landau gauge renormalization constants, [3U ESI USE EZl , and, [Hj IH1 ITu] 12*2*]. 
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in the MS scheme. This is useful since it checks that the expansion of the one loop diagrams 
has been performed correctly when the coupling constant and mass are replaced by their renor- 
malized variables. The quantity U' 1 (m 2 ,m 2 ) is itself clearly finite by simple power counting. 
Finally, we arrive at our expression for the two loop MS pole mass using the LCO Lagrangian 
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where 52 = (4v / 3/3)C/2(tt/3) and Cl2(x) is the Clausen function. 

Another check on the symbolic manipulation routines we have written was to consider the 
Curci- Ferrari Lagrangian, |19j . in the Landau gauge which is effectively QCD with a gluon mass 
included by hand. The corresponding Lagrangian is 
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where we have included the ghost mass term for completeness. As the one loop pole mass for 
this theory was given in |26l I25j . it does not require much more effort to produce the two loop 
MS correction using the same symbolic manipulation programmes. In this case the same Qgraf 
output is used but with a null a vertex and the usual quartic gluon interaction. The same check 
that the 2-point function is finite was satisfied. By contrast to 1)3. 7 Jl we find that the pole mass 
for (|3.8|) is, with s p 
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Although the Lagrangian for QCD in the non-linear Curci-Ferrari gauge formally differs from 
that for the Landau gauge in relation to the ghost gluon interaction term for a = 0, we have 
checked that the same pole mass emerges as (|3.9j) for the usual covariant Landau gauge fixing. 



4 Analysis. 



The next stage of our analysis is to produce an estimate for the pole mass. In [27| another 
method of estimating the Gross-Neveu mass gap was used. This required knowledge of the full 
2-point function as a function of the external momentum. This approach is not available to us 
for QCD since the technology does not exist to compute the gluon 2-point function exactly as a 
function of momentum. Instead we simply extend the argument of to two loops. In jS] an 
estimate for an effective gluon mass was obtained by examining the location of the minimum 
of V{a) using Q2.9|) . The effective gluon mass is essentially the bare mass of L a . It does not 
take account of quantum corrections. In t25 j it was argued that a more appropriate quantity 



to estimate from the effective potential was m^CO itself. Specifically the quantity V (ttv^qq) 
was constructed by inverting the one loop part of 1)3 ,7|) to obtain m 2 (//) as a function of m^QQ 
before substituting for m 2 (/i) using the relation with (a) 
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Using the minimization criterion, |25| . 
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emerged where we have not included the explicit expansion of £(<?). This is because it would 
introduce an unnecessary truncation error into the estimates for the pole mass. Ignoring the 
trivial solution of fn^QQ = which corresponds to the unstable vacuum, the non-trivial condition 
determines the pole mass estimate at one loop. To solve this the renormalization scale [i was 
parametrically related to "^CO ^Y m LCO = s ^ 2 which leaves a parametric relation between the 
running coupling constant and to^qq 



3QC A (16T F N f - 35C A ) [(3465tt\/3 + 46201n(s) - 25690 J C 2 A - 8UC F T F N f 

+ (19240 - 158471-^/3 - 3792 ln(s)) C A T F N f 



+ (7681n(a)-3328)T£jV| 



(4.5) 



where y = Ca5 ,2 /(167t 2 ). However, from the one loop /3-function the coupling constant can be 
related to the fundamental scale using 



where 



16vr 2 
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3 3 

Hence, the [i independent estimate for the pole mass emerged, |25j . 
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Equipped with the two loop pole mass of Q3.7JI . we now repeat the above analysis by first 
inverting mf^QQ to obtain m 2 (fj.) at two loops as a function of ^LCO' Substituting this into 
(|2.7j) and truncating at two loops, we find 
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where again we have not substituted for £(g). As this general expression is rather cumbersome, 
in order to illustrate the two loop analysis we concentrate for the moment on the case of SU(3) 
with Nf = when we simply have 
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(4.10) 



Solving 1)4.3)1 as before and discarding the trivial solution yields the two loop correction to 1)4.4)) 
for 577(3) Yang-Mills theory, 
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In analysing this along the lines of the one loop case, it transpires that the resulting two loop 
correction for the it^qq estimate is not fi independent. Therefore, we choose to return to the 
procedure of fSj and select the scale /j? so as to remove the logarithms in 1)4.11)1 . This fixes 
g{n) to a particular numerical value but using it the mass scale is recovered from the two loop 
extension of ()4.6)) 
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(4.13) 



N f 


1 loop 


2 loop 





2.10 


2.13 


2 


1.74 


2.21 


3 


1.55 


2.32 



Table 1. One and two loop estimates of mLCo/^Tru f° r SU(3). 



We have obtained estimates for both groups SU(2) and SU(3) for several quark flavours. 
These are summarized in Tables 1 and 2. Several features emerge. First, for Yang-Mills interest- 
ingly the two loop correction is less than 2% percent of the one loop value which suggests that 
our approximation is reliable. Unfortunately when quarks are included for both colour groups 
the situation is different with the two loop estimates being significantly larger than the one loop 
ones. Though for SU(3) they are of a similar size as the Yang-Mills value. Given the stability 
of the two loop results for Nt = compared with Nf ^ 0, it would suggest that the analysis 
when quarks are present is lacking some stabilising ingredient. One possibility is that for full 
QCD one actually requires quark masses to be included. 



13 



N f 


1 loop 


2 loop 





2.10 


2.13 


2 


1.54 


2.29 


3 


1.24 


2.58 



Table 2. One and two loop estimates of tolco /^WjiC for SU(2). 

5 Discussion. 

We have produced estimates for the gluon pole mass in QCD from the local composite operator 
method which systematically introduces an extra scalar field coupled to the gluon mass operator 
into the Lagrangian. The one loop renormalization scale independent estimate of [25] is stable 
to the two loop corrections for Yang-Mills theory but in the presence of quarks the estimates 
were significantly different. To improve the convergence for this case one could introduce masses 
for the quarks either by hand or by extending the LCO formalism to include the analogous mass 
operator vacuum expectation value {ipip) which is clearly beyond the scope of the present article. 
Moreover, if one accepts that a gluon mass emerges dynamically in QCD, one would then have 
to include gluon mass corrections in the estimates of the pole mass of the quarks. Although 
we have followed one procedure to deduce estimates for the gluon pole mass, other methods 
are possible. Indeed knowledge of the full momentum dependence of the gluon 2-point function 
would allow for the possibility of repeating the two loop analysis which was carried out in [22| 
using Grunberg's method of effective charges, [3*51 13*9*] . 
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